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' We present the results of complete tree level calculation for W boson 

. production processes e~p e~ ^^Vfj,X and e~p — > e' ^"v^X introduc- 
ing anomalous WWj and WW Z couplings. Detailed results for the 

Q^' distributions of final state particles are obtained. In the region of small 

ji^i momentum transfer we calculate the contribution of hadronlike photon 

Q ' component in the structure function approach. 



1 Introduction 

In recent years the charged and neutral current sectors of the Standard Model 
have been tested with excellent precision in the experiments at LEP and SLC 
However the gauge boson sector still remains practically untouched by direct 
measurements of high accuracy. Deviations of three and four gauge boson cou- 
plings from the Standard Model values would be an obvious signal of some new 
physics. 

At present time the best limits on anomalous three vector boson couplings 
are given by CDF and DO data (Fermilab Tevatron, a/s = 1.8 TeV) From the 
measurement of W'y and WW production these collaborations set the limits of 
order 1 on the deviations of (/c, A) couplings (see section 2.1) from the Standard 
Model values. Significant improvement of these limits (one order of magnitude) 
will be achieved by the detection of WW production at LEP2 0. 

In this paper we consider the possibilities of HERA ep collider (30 GeV elec- 
trons on 820 GeV protons, ^/s = 314 GeV) for the measurement of the vector 
boson anomalous couplings. The main difference between our study and the 
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previous investigations is the exact calculation of tree level amplitude for com- 
plete set of diagrams with the four particle final state, including nonstandard 
WW'j and WWZ vertices. At present time the luminosity of HERA (several 
pb~^ /year) is too small to produce sufficient number of W bosons. However, 
after the luminosity upgrade to 100 - 200 pb~^ the detection of anomalous signal 
or setting new limits on the anomalous couplings becomes realistic. 

2 The reactions e~p e~/i+z/^X, e~jj.~u^X 

It is known from the previous study (in particular we would like to distinguish 
the paper [^) that in the Standard Model W^-bosons are produced in ep scatter- 
ing mainly in the channels ep —>■ eWX. The contribution of the channel with 
neutrino in the final state ep — >■ V(,W~X is 20 times smaller. The following de- 
cay of W boson to muon and muonic neutrino produces the four-fermion state 
e~/iz/^g, and the corresponding event signature is muon(antimuon) with missing 
transverse momentum. The signal of W boson production in leptonic channels 
can be observed easier than in hadronic W decay channels, where large QCD 
background processes must be carefully separated from the signal. 

Ten Feynman diagrams for the reaction e~qi — >■ e~/i^z/^g2 are shown in Fig.l. 
All diagrams have intermediate W and the properties of final state are defined 
by W interaction dynamics. In this sense all diagrams are W+ producing' and 
there are no irreducible background graphs that could be neglected in order 
to simplify the procedure. If we replace diagram 5 by similar one where W~ 
boson is radiated from the initial electron and change /x^z/^ to /i~z/^, gi,2 ^1,2, 
we obtain a set of ten diagrams for the process e~p — > e^fi^u^X which is 
W~ producing'. Diagrams 4 and 10 in both cases contain WW'y and WWZ 
vertices and in the following we shall use for them nonstandard gauge invariant 
structure. 

If we separate subsets of diagrams from the complete tree level set in Fig.l 
and then separate Feynman subgraphs from these subsets, we obtain some ap- 
proximations that were used in the previous calculations of W production pro- 
cesses. The simplest approximation is given by diagrams 2,3,4 where t-channel 
photon and s-channel W are taken on-shell (7*^1 — >■ W*q2, 2 — > 2 subpro- 
cess approximation). If we integrate then with equivalent photon structure 
function for incoming 7* and consider the decay of on-shell W to /xi/^, rather 
satisfactory estimate of total cross section can be obtained. The calculation in 
7*gi — * W*q2 subprocess approximation with anomalous vector boson couplings 
can be found in Weak correction to this result is given by diagrams 8,9,10 
(containing subprocess Z*qi — >• W~^q2). In the case if the accuracy of equiva- 
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lent photon approximation is not sufficient, at the next step one could consider 
e~qi e~W^q2, i.e. 2 — 3 process approximation with on-shell vector boson. 
Complete tree level calculation for the process e^qi e~W~^q2 by means of 
helicity amplitude method was performed in [Q. If we take the amplitude with 
W off-shell decaying to fermionic pair (for instance, e~qi e~^~^i'^q2 in the 
2 — >• 4 process approximation), the subset of diagrams (2,3,4,5,8,9,10) becomes 
gauge noninvariant and in order to restore the gauge invariance it is necessary 
to add ladder diagrams 1,6,7. In the papers || complete tree level calculation 
was performed for the case of Standard Model. We generalize these calculations 
in the case of complete tree level 2 — > 4 muonic channels with anomalous C and 
P conserving three vector boson couplings, and compare some of our results 
with 10] where similar analysis was done by means of EP VEC generator . 
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Figure 1: Feynman diagrams for the process e qi 



While the results for total rate provided by 2 — 2 and 2 — > 3 approxima- 
tions can be quite satisfactory, it is not possible to calculate precisely most of 
the distributions of experimental interest. The accuracy of equivalent photon 
approximation becomes rather poor (see section 4) especially at large transverse 
momenta and the narrow-width approximation for the W is usually not good 
near the W production threshold. At the same time it is obviously difficult to 
make any conclusions about the origin of new phenomena observing only the 
deviation of event counting rate from the Standard Model value. It is important 
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to know what regions of phase space are affected by new interaction dynamics 
and what is the ratio new signal/ background, i.e. to calculate precisely the 
distributions of particles in the final state. 

Largest contribution to the cross section of e~qi e~/i+t'^g2 process is 
given by diagram 3 when photon and t-channel quark are close to mass shell 
In this configuration QCD corrections become large and t-channel interme- 
diate quark can appear nonperturbatively as a constituent of photon ('resolved 
photon' contribution) when it is usually described by experimentally measured 
gamma structure function. The process of W production in this picture is quark- 
antiquark q^qp fusion to W when photon fragments into quark constituents q^ 
before interacting with the proton constituent quark qp. It was shown in 
that resolved photon mechanism is not dominant, but nevertheless the region 
of small t requires special consideration and careful separation of 'resolved' and 
ordinary contributions is needed. To be sure that the numbers are not changed 
significantly by the new parametrizations of 7 and p structure functions, we re- 
peat the cross section calculation of the resolved part using the scheme similar 
to one proposed in [Q. 

2.1 Anomalous three vector boson couplings 

General effective lagrangian of two charged and one neutral gauge boson in- 
teraction was proposed in 0. The restrictions on the lagrangian imposed by 
invariance under discrete symmetries and gauge invariance were considered in 
[ p!0[] . U(l) gauge invariant, C and P parity conserving effective lagrangian has 
the form 



Leff = gviW+W^'V - W^^'W+V, + k W+W.V^" + -^w+w^,v^) (1) 

where = e and gz = ecos^w / sin^w , = d^W^ - d^Wf,, V^y = d^V^ - 
duVfj_. Spatial structure of the fourth term in the lagrangian of dimension six 
is multiplied by factor, so parameters A and k are dimensionless. In the 
momentum space if all momenta are incoming {pi + P2 + Pa = 0), we have the 
following expression for W~^{pi)W~{p2)V{p3) vertex 



^f,upiPi,P2,P3) = 9v[9t,v{Pi -P2 —{{P2Pz)Pi - {pm)P2)p (2) 

+9pp{kpz -Pi + -^{{P2P3)pi - {PlP2)Pz)u 
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+9up{p2 - kp3 ^((P1P3)P2 - (PiP2)P3)m 

+ ^iP2uP3uPlp - P3uPluP2p)] 



m 



w 



In the special case X = 0, k = 1 this vertex reduces to Standard Model one. 



3 Resolved photon contribution 

In this section we follow the scheme proposed in 0] for the calculation of resolved 
photon contribution (diagram 3, Fig.l), but in our case W boson is off-shell and 
consequently the scale of the equivalent photon approximation is different. Sim- 
ilar procedure (including the corrections from 'finite terms') was considered in 
[0]. We separate the 'resolved' and 'direct' production mechanisms at the scale 
= (gi — qwY ill and qw are the momenta of initial quark and intermediate 
W, correspondingly): 

O' — (Jresolved + J ~~d^ 

Second term in this formula ('direct') is calculated numerically using the exact 
matrix element for 10 diagrams (see Fig.l) with the transferred momentum 
cutoff at t = A^. Resolved photon cross section in the case of monoenergetic 
initial gamma on shell is given by the convolution of photon structure function 
(measured in 77 collisions) and q-^Qp ^ W fusion cross section 

O-resolved = J dxidx2fqi/'y{xi, Q^) a{qiq2 ^ W) fq^/p{x2, Q^) (4) 

and Qp are the scales for photon and proton structure functions, correspond- 
ingly. Using the Breit-Wigner formula in the approximation of infinitely small 
W width we get 

cr(gig2 ^W) = ^^^^GFm^\Vi2\^6{xiX2S - m^) (5) 

where V12 is the CKM matrix element for charged current ipq^ipq^] here and 
in the following we are not indicating the sum over possible quark spieces. 
The experimentally measured photon structure function fq^/-f includes point- 
like as well as hadron-like parts. In the leading logarithmic approximation the 
perturbative (point-like) part of photon structure function can be expressed as 

f^lM^ Q?) = ^[^' + (1 - (6) 
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where Cg is the quark charge. This part was aheady taken into account by our 
calculation for direct contribution in (3). In order to avoid double counting in 
the contributions from gamma structure function and from direct process it is 
necessary to subtract the point-like term (6) from /^^/^(xi, Q"^). This procedure 
was illustrated explicitly for the more simple example of the reaction 7^1 Vq2 
{y = Z, W) in [0 and it was shown that in the case of Z production indeed the 
LO counterterm (6) rescales the leading logarithmic structure in the direct part 
from to Q'^ when /^^/^(xi, Q^) is taken at momentum transfer Q'^. Finally the 
dependence from is absent in the sum of resolved and direct contributions. 
Our case is certainly more complicated but double counting can be avoided at 
least on the leading logarithmic level. 

Introducing the usual equivalent photon approximation |TT| for gamma in 
the initial state 

fq/e{XjQ^^Qww) ~ [ f gi/'y{~ J Qy) fy/eiy, Qww) (7) 
J X V V 



where 



dMv, Q'ww) = |^[^^^^ - 2^e^]^Q' (8) 



Qmin = mlx'^/^l — x) and Q^ax = Qww defined by some process scale 
that will be discussed later. After the substitution of (5) and (7) into (4) and 
subtraction of counterterm (6) we finally get 



f-y/eiV) Qww)fq2/pi ^ ^ 5 Qp) 

There are four scales QiY]y,Q'^,A'^,Qp to be defined in this formula. The 
Weizsacker- Williams scale Q^^y can be chosen equal to A^. It is easy to jus- 
tify this choice [§] looking at the distribution da / dlogQ^^^/ (Fig-2) calculated 
exactly for ten diagrams of 2 — 4 process in Fig.l. The flat part of this dis- 
tribution corresponds to the cross section behaviour da/dt ~ 1/t (not as 
it seems at the first sight, double poles are cancelled), and rapidly decreases 
starting from A^. We always take proton structure function scale = m^. 
The values of Q^,A^ are arbitrary and the final result for the sum of resolved 
and direct contributions (3) should not depend essentially from the choice of 
these two scales. 

In Table 1 we show the results for resolved photon cross section (9) calcu- 



lated by means of BASES MC integrator |12], using proton structure functions 
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MRS A 


CTEQ3m 


Qww 




A 


DGl 


LAC2 


GRV LO 


DGl 


LAC2 


GRV LO 


0.2 


mw 


0.2 


-11.6 


-4.1 


-7.3 


-11.9 


-4.2 


-7.6 


0.2 


m^y/lO 


0.2 


-7.5 


2.0 


-3.3 


-7.6 


2.7 


-3.1 


1.0 


mw 


1.0 


-5.6 


3.3 


-0.6 


-5.7 


4.0 


-0.1 


1.0 




1.0 


1.7 


13.4 


6.5 


1.7 


14.8 


7.4 


5.0 




5.0 


4.3 


14.7 


9.8 


4.4 


16.3 


11.0 


5.0 


mw/W 


5.0 


15.7 


29.6 


21.2 


16.0 


32.0 


23.1 



Table 1: Resolved photon cross section of the process e~p e~ jj.'^v^X in fb 
(see formula (9)). Different sets of photon and proton structure functions were 
used. 



MRS [13 1 and CTEQ |14 



151, LAC-G 16 



and GRV [17 



photon structure functions DG 
Last two versions of gamma structure function are improved 
parametrizations in the framework of the approach [p!5| ]. 

One can see that if proton structure function is measured with rather good 
accuracy and two parametrizations we are using give similar results (consistent 
within the one standard deviation error of our MC integration), photon struc- 
ture function is still poorly known. LAC-G parametrization gives the cross 
section regularly smaller than the values obtained with the help of DG and 
GRV parametrizations. We checked that our point obtained with the help of 
LAC2 and HMRS B structure functions at = m^/10, A = 5 GeV and equal 
to aresoived = 22 fb, is close to the value 24 fb given in ^ at the same parameter 
values. So the possible correction to our result (9) from the so-called 'finite 
terms' is around 8%, which is much less than the difference of results obtained 
by using different photon structure functons. Negative cross section in Table 1 
means that at a given scale of A strong double counting regime of direct and 
resolved contrubutions takes place. In other words, most part of the resolved 
cross section is already taken into account by the calculation for direct term in 
(3). 
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4 Complete tree level calculation and anoma- 
lous signal of W in the distribitions 

4.1 General framework 

Complete tree level calculation of 10 diagrams in Fig.l with anomalous three 
vector boson couplings (direct process) and the following generation of particle 



distributions were done by means of CompHEP package [|l^, |T9|. The ampli- 
tude corresponding to 55 squared diagrams and interferences between diagrams 
was calculated symbolically. In order to avoid t-channel poles, masses of elec- 
tron and quarks were kept nonzero. [| Equivalent photon approximation was 
not used. After that symbolic expressions are automatically converted to FOR- 
TRAN codes and linked to special program for seven dimensional Monte-Carlo 
integration over four particle phase space, and adaptive integration package VE- 



GAS |]21|. In the process of four particle phase space generation we introduce 



so-called kinematical regularization of the peaks inherent to the amplitude 
under consideration. Especially these are t-channel gamma peaks in the ee'j 
vertices of diagrams 1-4 (Fig.l), ^fi'j vertex of diagram 1 and py-resonance 
peak in the diagrams with s-channel VT-boson. 

CompHEP package is a software product in the framework of one of a few 
general approaches [18, 2^ developed in recent time for the analysis of multi- 



particle exclusive states at new colliders, when hundreds of Feynman diagrams 
contribute to the amplitude and should be exactly calculated. More details can 
be found in |]T^ . 



In our calculations we used the Breit-Wigner propagator with constant width 
for the ly-boson. Generally speaking, if we have some complete tree level set of 
diagrams, straightforward replacement of lowest order vector boson propagator 
by the propagator with finite width violates gauge invariance of the amplitude 
and can break gauge cancellations between diagrams, leading to numerically 
unstable false results. For this reason we used the well-known "overall" form of 
propagator replacement when the entire amplitude is multiplied by a factor 

2 2 



Generally speaking in other cases this prescription could affect strongly nonres- 
onant terms in the amplitude p3[ (in the region of phase space where ~ ^w)-, 



^YoT this reason we need A cut only for matching of direct and resolved parts of cross 
section. For instance, if m„ — 5 MeV, rud = 10 MeV, nis = 0.2 GeV and TOc =1.3 GeV and 
there are no kinematical cuts, adiriep — > e/i"'"t'^X) = 102.5(5) fb in the standard case with 
MRS A structure functions. 
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but the case under consideration is free from this difficulty. 



A = 0.2 GeV 


A 


k 


eu — > efi^Ufj^d 


ed efi^u^u 


eu — > efi'^u^s 


es efi^u^c 


(^tot 





1 


68.2(4) 


15.2(1) 


3.5(0) 


5.2(0) 


92.1(5) 


1 


1 


72.6(5) 


15.6(1) 


3.7(0) 


5.3(0) 


97.2(6) 








52.6(5) 


13.6(1) 


2.7(0) 


4.6(0) 


73.5(6) 





2 


99.6(6) 


17.2(2) 


5.1(0) 


6.3(0) 


128.2(8) 


A = 


L.O GeV 





1 


61.2(4) 


12.4(2) 


3.1(0) 


5.2(0) 


81.9(6) 


1 


1 


66.0(4) 


12.6(1) 


3.4(0) 


5.3(0) 


87.3(5) 








45.3(3) 


10.8(1) 


2.3(0) 


4.6(0) 


63.0(4) 





2 


92.7(6) 


14.9(1) 


4.8(0) 


6.3(0) 


118.7(7) 


A = 5.0 GeV 





1 


52.5(3) 


8.9(0) 


2.7(0) 


4.5(0) 


68.6(3) 


1 


1 


57.4(4) 


9.0(0) 


2.9(0) 


4.6(0) 


73.9(4) 








36.4(3) 


7.4(0) 


1.9(0) 


3.9(0) 


49.6(3) 





2 


83.8(5) 


11.1(0) 


4.3(0) 


5.6(0) 


104.8(5) 



Table 2: Total cross sections (fb) of the main partonic producing processes 
(see 10 diagrams in Fig.l) in the reaction e~p — > e~ ^'^v^X for the various sets 
of /c, A. Invariant kinematical cut A^ = [pq —pwY- Proton structure function 
MRS A. One standard deviation error of Monte-Carlo integration for the last 
digit is indicated in brackets. 



The accuracy of our Monte Carlo calculation of adir is usually 0.6-0.8% (see 
Table 2). This choice is related to the precision of proton structure function 
parametrization. We checked that if we replace the proton structure set MRS A 
13 that we are usmg by the set CTEQ3m IH, the relat ve difference of results 



obtained with these two sets does not exceed 1%. 

Four partonic processes from 12 possible give the main contribution to W 
boson production. We present the results of total cross section calculation for 
the main subprocesses of the production channel e~p — > e~ ii^v^X in Ta- 
ble 2. Main contribution comes from the subprocess e~u e^fi^u^d. The 
remaining 8 partonic reactions have very small individual total cross sections 
and the sum of their contributions is of order 1 fb. The cross section of W~ 
production channel e~p — >■ e'n'u^X is compared with the case of produc- 
tion in Table 4. In the case of W~ production the main partonic subprocess is 
e~d — > e~fj,~^h'^u and the total rate is slightly smaller because there are less d- 



9 



A = 0.2 GeV, Ef, > 10 GeV, missing pr > 20 GeV 


A 


k 


eu — > efi^v^d 


ed — > eji^'u^j^u 


eu — > e/jL'^i'fj^s 


es — > eji^v^c 





1 




1 
1 

2 


49.2(4) 
54.7(5) 
36.6(4) 
75.9(5) 


11.1(1) 
11.4(1) 
9.8(1) 
12.8(2) 


2.5(0) 
2.8(0) 
1.8(0) 
3.9(0) 


3.6(0) 
3.7(0) 
3.0(0) 
4.5(0) 


66.4(5) 
72.6(6) 
51.2(5) 
97.1(7) 



Table 3: The same as in Table 2 with the kinematical cuts imposed > 10 
GeV, missing > 20 GeV. These cuts are used to exclude the misidentification 
backgounds in the electron channel (see section 4.2). 

quarks in the proton than w-quarks. Resolved photon contribution from Table 1 
is also indicated. One can see that indeed rather weak dependence of adir + Cres 
from the cutoff A takes place in so far as the decrease of direct part with the 
growth of A is compensated by the increase of resolved part. 

We already mentioned in Section 3 that — was always taken as the 
momentum transferred scale in the proton structure functions. It is important 
to find out how the changes of affect the total cross section value and compare 
the possible deviation of total rate caused by the change of hadronic scale with 
the deviation of total rate coming from anomalous three vector boson couplings. 
The uncertainties coming from the value of should be less than the effect of 
anomalous couplings to make the phenomenological restrictions based on the 
value of total rate more meaningful. We show the numbers for total cross section 
of the main producing partonic process eu — > efiv^d calculated for proton 
structure function scales m^/2 and 2m^ in Table 5. These values are taken 
as illustrative ones because in real partonic processes the contributions from so 
small/large values of are negligible. We can see that the total cross section 
deviation from the standard choice of hadronic scale Q"^ = is around 2.5 - 
3%, while the effects of anomalous W couplings are much larger. For instance, 
the effect coming from anomalous A;-term in 20 - 30% and the effect of anomalous 
A term is 6 - 8% (Table 2). 

4.2 The reactions e~p — ^ e'e^VeX, e~p — ^ e~e~UeX and the 
comparison with EPVEC generator 

Total number of diagrams in the channels e~p e'e'^v^X and e~p e~e''DeX 
is 20 (ten additional diagrams to Fig.l with i-channel 1^-boson exchanges ap- 
pear). However at HERA energy the correction coming from these additional 
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92.1 


81.9 


68.6 
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1 


97.2 


87.3 


73.9 










73.5 


63.0 


49.6 
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128.2 


118.7 


104.8 






-4.1 


3.3 


14.7 


O-dir + (^res, W+ , SM 




88.0 


85.2 


83.3 







1 


80.3 


68.6 


52.4 


TdiriW-) 


1 


1 


82.7 


70.6 


53.9 










69.5 


57.6 


41.4 







2 


97.6 


86.0 


69.5 






-11.0 


-3.8 


7.9 


CTdir + (^res, W' , SM 




69.3 


64.8 


60.3 



Table 4: Total cross sections (fb) of the reactions e~p e^fi^u^X {W^ pro- 
duction) and e~p — >• e~/i~z/^X (W~ production) in the case of standard (SM) 
and anomalous three vector boson interaction at different values of A cutoff 
parameter. Proton structure function MRS A. In the calculation of resolved 
photon contribution gamma structure function LAC-G2 was taken at the scale 
mw ■ 

diagrams is very small because they do not contain t-channel photon poles. 
Usually this weak contribution is neglected. Such approximation was used in 
the recent simulation for HERA (by means of EPVEC generator f^), where 
event topology and realistic kinematical cuts were considered in more details. 

The electron channels e'e^UgX and e~e~z/eX are more difficult for exper- 
imental study than the muon channels e'fi^u^X and e'fi^u^X. First, large 
background from the neutral current deep inelastic scattering process ep eX 
appears in the case when the final jet energy is not completely registered in the 
hadronic calorimeter and for this reason some missing pt is observed. Second, 
large background from the charged current DIS process ep —>■ vX appears in the 
case when 7r° from the final jet is misidentified as e~{e^) giving again the final 
state with e~(e"*"), jet and missing px- In order to suppress these misidentifica- 
tion backgrounds the following kinematical cuts were used in (1) isolated 
electron with the energy Ee > 10 GeV (2) missing transverse momentum px > 
20 GeV are required. The requirement of isolated electromagnetic cluster is 
removed for the case of muonic channels, and 10 GeV energy cut seems too 
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strong for the final muon, since good reconstruction of several GeV muons is 
available in ZEUS and HI detectors. Missing transverse momentum cut also 
seems not necessary for the muon channels (scattered lepton is different from 
W decay lepton). 

We checked by direct calculation that in the absence of missing px cut, muon 
energy cut at 2-3 GeV practically does not affect the numbers for total cross 
sections shown in Table 2 and Table 4. The combination of muon energy cut at 
10 GeV and missing pt cut at 20 GeV gives the cross sections by 30-35% smaller 
than before cuts (see Table 3). Following [0 we shall consider electron channels 
in the approximation of 10 diagrams subset (Fig.l) with the kinematical cuts 
Ee > 10 GeV and missing pr > 20 GeV. 

In Fig. 3 we show the distributions of energy, angle with the beam and trans- 
verse momentum for electron, muon and final quark in the muonic channel 
e^^^u^X, obtained by means of CompHEP (Standard Model case). The same 
distributions for electron channel, obtained with the help of EPVEC generator, 
can be found in 0. (We are using in (3) A = 0.2 GeV, while the distribu- 
tions in 1^ are calculated using A = 5 GeV.) Although the normalization is 
not indicated in [0, one can observe that the agreement of the shapes is sat- 
isfactory. Soft muons in the distributions da/dE^ and da/dpT^ come from the 
ladder diagrams 1,6,7 in Fig.l. Jets at the angle 180 degrees with the proton 
beam appear from diagram 3, Fig.l, when the quasireal photon produces two 
quarks coUinear to initial electron. Soft muons and backward jets are absent in 
[0], because besides the E^ and missing pt cuts mentioned above, the EPVEC 
generator contains build-in cuts [^] separating some region of phase space near 
the W pole. 

4.3 Sensitivity to anomalous couplings 

Let us return to our calculation with anomalous interaction of vector bosons. 
We can estimate approximately the possibilities of HERA in the detection of 



anomalous couplings using a simple criteria (see, for instance, [^) for the 
number of events that is necessary to observe Act deviation from the total 
cross section value a: 

(11) 

It follows from Tables 2,4 that at the integrated luminosity of HERA L = 
200 pb~^ in the channels of W'^ and W~ production e~ ii'^v^q and e~ ii~u^q 
we shall have about 35 events/year. Deviation of A in the vertex (1) from the 
zero standard value AA = 1 gives us 5% deviation in the total W production 
rate and we need about 400 events to observe it. However the deviation AA;=1 
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(AA; = A; — 1) changes the total cross section by 20-40 % and less than 25 events 
will be needed for some experimental evidence. and W~ production in 

the channels e'e'^u^q and e~e~Ueq will give slightly less reliably reconstructed 
events than in muonic channels (kinematical cuts must be introduced to tune 
off the misidentification backgrounds), so in total around 60 W bosons/year 
decaying to electrons and muons could be observed. It follows that it will be 
difficult to improve CDF and DO limits on A [Q, but some improvenent of Ak 
restriction could be possible. 

More precisely, systematic errors on the detector acceptance A and the un- 
certainty in the luminosity measurement L should be taken into account. The 
former are estimated on the level of 2%(1%) for the integrated luminosity of 
order 10^ pb~^ {10^ pb~^) and the latter is taken to be 2%. The uncertainty of 
the total cross section measurement has the form ^ 

Act , 1 ,AL,, ,AA,^, i , , 

The acceptances in both electron and muon channels are taken to be 65%. From 
equation (12) we derive the following limits for Ak and A, giving the observable 
deviation of total cross section from the Standard Model value at 68% and 95% 
confidence level: 

-1.70 < A < 1.70, - 1.05 < AA; < 0.48, 68% CL 
-2.24 < A < 2.24, AA; < 0.89, 95% CL 

at the integrated luminosity 200 pb~^, and 

-1.03 < A < 1.03, - 0.31 < AA; < 0.27, 68% CL 
-1.75 < A < 1.75, - 0.58 < AA; < 0.46, 95% CL 

at the integrated luminosity 1000 pb^^. Here only one coupling from the pair 
(A, A;) is assumed to be different from the SM value. We do not indicate the 
negative 95%CL AA; limit at the integrated luminosity 200 pb~^ because the 
cross section deviation from the SM value stops to increase starting from AA; ~ 
-1.5 and the effect cannot be observed with small statistics. However, positive 
AA; limit at the same luminosity is competitive with the early expectations from 
LEP2 0. These limits could be of course improved by taking into account the 
channels W — > jets with final electron and three jets in the final state. Low 
acceptance in the jets channel (~ 20%, |^) and complicated situation with the 
separation of QCD backgrounds requires an independent careful study, and we 
are not considering this possibility here. 
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Of course the calculation of total rate is very important and the ratio 
CtotiW) I atot{Z) that was considered in |Q could be the clear indicator to anoma- 
lous gauge boson coupling. However the only way to see definetly if the devi- 
ation of the ratio is really due to anomalous WW'j interaction but not caused 
by some other reason, is to inspect what regions of phase space are affected by 
anomalous W interaction dynamics and how they are affected. It is natural to 
use, as proposed in [H, the distributions of final quark and muon transverse 
momenta. 

We show the distributions of final quark jet transverse momentum in Fig. 4, 5. 
It follows from symbolic calculation in 2 — 2 approximation that the cross 
section depends from A quadratically 0]. We checked at complete tree level 
2 — >• 4 that this is true for atot at about 1% accuracy and no difference at 
positive and negative values of A is observed in the distributions. The deviation 
of k is clearly seen (Fig. 5), dcr/dpT becomes harder when k is less than standard 
value k = 1. 

The distributions of final muon transverse momentum are shown in Fig. 6, 7. 
Similar to the previous case of jet pr distribution, the effect coming from A 
is very small and the dependence of distribution from k is rather strong. The 
forward and backward slopes of W peak can be slightly shifted if we take into 
account W production by the resolved photon (see the details inQ), but this 
shift is the same for standard and nonstandard cases. 

The distributions of final muon rapidity for the standard and anomalous 
cases are shown in Fig.8. The direction of proton beam was chosen as the 
direction of positive rapidity axis. 

An important point concerning the Pt distributions of the quark and muon 
for the direct process is their sensitivity to higher order QCD corrections, which 
could be large. As it was stated in [|, ||, the integration close to w-channel pole 
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1 


70.6(5) 


68.2(4) 


66.3(6) 
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75.5(5) 


72.6(5) 


71.1(4) 








54.1(4) 


52.6(5) 


51.2(4) 





2 


102.3(6) 


99.6(6) 


96.2(6) 



Table 5: Total cross section (fb) of the main producing partonic process 
eu — i> efiUf^d calculated using three scales in the parametrizations of proton 
structure function MRS A. 
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A 


k 


exact result 


QwwA — fnw 


QwwA='2,^ GeV 





1 


61.2(4) 


66.6(8) 


61.1(5) 


1 


1 


66.0(4) 


70.9(6) 


64.5(5) 








45.3(3) 


49.4(5) 


46.5(4) 





2 


92.7(6) 


96.7(9) 


88.1(8) 



Table 6: Comparison of exact calculation (fb) for the process e~M — e~ ^'^v^d 
with equivalent photon approximation for the process 7M ^Ji^v^d calculated 
at the scale Qww (8) equal to niw and 20 GeV, A =1.0 GeV. 

in the diagram 3, Fig.l involves the momenta of the order Aqcd in the small 
Pt region, when the QCD corrections can be expected to be significant. Total 
cross section in the absence of w-channel or cuts contains some degree of 
uncertainty. It was mentioned in that in connection with normalization 
uncertainty of the cross section the events when the jet escapes detection and 
only the lepton with missing pt are observed could provide an important test 
on the normalisation of W production rate. 

Finally we would like to discuss the question of equivalent photon approx- 
imation ||ll| accuracy in our case. In Table 6 we compare exact result for the 
partonic subprocess eu efi'^u^d and the equivalent gamma approximation 
calculation for the process 7M — ^'^v^d. In the latter case intial photon mo- 
mentum is distributed according to (8) where we used two different choices of 
Qww- One can see that the 'natural' choice Qww = ^w overestimates the 
cross section by 8% (let us remind (see Table 4) that the effect of AA = 1 is 
5%). It is possible of course to adjust Qww which is not strictly fixed at any 
value, but defined by some typical process dependent momentum transferred 
scale, and get agreement of exact and WW cross sections for the standard case 
A = 0, /c = 1; it turns out that the corresponding value is Qww = 20 GeV. 
However, after fixing of this value in the anomalous case X = 0, k = 2 again we 
observe 5% deviation. Equivalent approximaions become too rough if precise 
separation of the signal is needed. 

5 Conclusion 

We presented the results of complete tree level calculation for the W boson 
production processes at the energy of HERA collider, introducing anomalous 
C and P conserving three vector boson couplings (1). The main and W~ 
production channels e~p e~Z"''z//X and e~p e~l~i)iX {I = e, /x)were con- 
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sidered. Following the earlier publications p, Q we separated the phase space 
at some scale of momentum transferred from the constituent quark to VT-boson 
in order to take into account the resolved photon contribution to the total rate. 
Resolved photon part was calculated in the structure function approach, using 
new parametrizations of photon and proton distribution functions. Perturba- 
tive (direct) part of the cross section was considered by means of CompHEP 
package |jl8[,lT^, when the tree level 2 —> 4 amplitude, corresponding to ten 



diagrams for each of and W~ production processes, is calculated exactly 
without any approximations. Some uncertainty in the normalization of total 
cross section exists for the reason of possibly large QCD corrections in the 
phase space region near the w-channel quark pole. In the muonic channels un- 
der consideration the total cross section is equal approximately to 150-160 fb, 
giving about 35 events/year at the integrated luminosity 200 pb~^. Kinematical 
cuts are necessary in the electron channels for separation of misidentification 
backgrounds, and the number of identifiable events from W eug is slightly 
smaller. 

We show explicitly what regions of phase space are affected by anomalous 
three vector boson interaction dynamics. In particular it follows from our analy- 
sis that even at the integrated luminosity 1000 pb~^ it will be extremely difficult 
to separate anomalous A term effect in (1) when A A is less than 1.5 (as already 
restricted by Tevatron data), but rather easy to observe anomalous k term ef- 
fect, when Ak is of order 0.4-0.8, which is strongly competitive with LEP2 
possibilities. 
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Figure 2: Distribution over logarithm of momentum transferred squared Qis = 
ti3 = (pg" — Flat part of the distribution corresponds to the behaviour 

da/dt ~ l/ti3 
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Figure 3: First row of plots - distributions of the electron energy, scattering an- 
gle and transverse momentum in the process e~p — * e~ ^'^v^X. Second row of 
plots - distributions of the muon energy, muon scattering angle and transverse 
momentum. Third row of plots - distributions of the quark energy, angle and 
transverse momentum for the same process. No kinematical cuts, all calcula- 
tions were done by means of CompHEP [18] (A = 0.2 GeV) to be compared 
with the same distributions obtained by means of EPVEC generator (A = 5.0 
GeV), see [7] 20 
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Figure 4: Distribution of jet transverse momentum in the reaction e~p — > 
e~ ix^v^X. Upper plot: no kinematical cuts, solid lines - standard case, dash 
lines - anomalous three vector boson couplings case, A = 1, A; = 1 Lower plot: 
the same distributioiis after kinematical cuts > 10 GeV, missing pt > 20 
GeV. 
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Figure 5: Distribution of jet transverse momentum in the reaction e~p 
e~fi~^i'^X. Upper plot: no kinematical cuts, solid lines - standard case, dash 
lines - anomalous three vector boson couplings case, X = 0, k = 0, thin solid 
lines - A = 0, k — 2. Lower plot: the same distributions after kinematical cuts 
> 10 GeV, missing pr > 20 GeV. 
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Figure 6: Distribution of muon transverse momentum in the reaction e~p 
e~ ii^v^X. Upper plot: no kinematical cuts, solid lines - standard case, dash 
lines - anomalous three vector boson couplings case, A = 1, A; = 1 Lower plot: 
the same distributioiis after kinematical cuts > 10 GeV, missing pt > 20 
GeV. 
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Figure 7: Distribution of muon transverse momentum in the reaction e~p 
e^fi^u^X. Upper plot: no kinematical cuts, solid lines - standard case, dash 
lines - anomalous three vector boson couplings case, X = 0, k = 0, thin solid 
lines - A = 0, k — 2. Lower plot: the same distributions after kinematical cuts 
> 10 GeV, missing pr > 20 GeV. 
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Figure 8: Distribution of muon rapidity in the reaction e~p — > e~ ji^ffj^X. Upper 
plot: no kinematical cuts, solid lines - standard case, dash lines - anomalous 
three vector boson couplings case, A = 1, A; = 1. Lower plot: no kinematical 
cuts, sohd lines - standard case, dash hnes - anomalous three vector boson 
couphngs case, A = 0, A; = 0, thin solid hnes - A = 0, A; = 2. 
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